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ABSTRACT 



Context. The surface of the Sun provides us with a unique and very detailed view of turbulent stellar convection. Studying its 
dynamics can therefore help us make significant progress in stellar convection modelling. Many features of solar surface turbulence 
like the supergranulation are still poorly understood. 

Aims. The aim of this work is to give new observational constraints on these flows by determining the horizontal scale dependence 
of the velocity and intensity fields, as represented by their power spectra, and to offer some theoretical guidelines to interpret these 
spectra. 

Methods. We use long time-series of images taken by the Solar Optical Telescope (SOT) on board the Hinode satellite; we reconstruct 
both horizontal (by granule tracking) and vertical (by Doppler effect) velocity fields in a field-of-view of ~ 75x75 Mm 2 . The dynamics 
in the subgranulation range can be investigated with unprecedented precision thanks to the absence of seeing effects and the use of 
the modulation transfer function of SOT for correcting the spectra. 

Results. At small subgranulation scales down to 0.4 Mm the spectral density of kinetic energy associated with vertical motions 
exhibits a £~ 10/3 -like power law, while the intensity fluctuation spectrum follows either a k~ [1/3 or a £~ 3 -like power law at the two 
continuum levels investigated (525 and 450 nm respectively). We discuss the possible physical origin of these scalings and interpret 
the combined presence of k~ n,i and k~ w/3 power laws for the intensity and vertical velocity as a signature of buoyancy-driven 
turbulent dynamics in a strongly thermally diffusive regime. In the mesogranulation range and up to a scale of 25 Mm, we find that 
the amplitude of the vertical velocity field decreases like A~ 3/2 with the horizontal scale A. This behaviour corresponds to a k 2 spectral 
power law. Still in the 2.5-10 Mm mesoscale range, we find that intensity fluctuations in the blue continuum also follow a k 2 power 
law. In passing we show that granule tracking cannot sample scales below 2.5 Mm. We finally further confirm the presence of a 
significant supergranulation energy peak at 30 Mm in the horizontal velocity power spectrum and show that the emergence of a pore 
erases this spectral peak. We tentatively estimate the scale height of the vertical velocity field in the supergranulation range and find 
1 Mm; this value suggests that supergranulation flows are shallow. 

Key words. Convection - Turbulence - Sun: photosphere 



1. Introduction 

A complete understanding of thermal convection in stars remains 
one of the main challenges of present day astrophysics. The 
Sun offers an unsurpassed detailed view to address this ques- 
tion. Solar convection is highly non-linear - typical Reynolds 
numbers are over 10 10 - making fully detailed direct numeri- 
cal simulations down to the viscous dissipation scales unafford- 
able. In spite of this difficulty, some observed features of so- 
lar convection like for instance granul ation are rather well un- 
derstood (e.g. Stei n & Nordlundl uOOO). Dynamical features at 
larger scales like the supergranulation pattern are much less un- 
derstood. Supergranulation is characterised by a horizontal ve- 
locity field spanning scal es from 15 Mm to 75 Mm accord- 
ing to the recent work of iRieutord et al .] (120081) . Despite sev- 
eral attempts, supergranulation has not be en identified in large- 
eddy or direct num e rical simulations (e.g. IRieutord et al.| [2002: 
Rincon et al. 2005; IStein et al.l |2009). Its origin remains un- 
known, althou gh some scenarios have been proposed. The clas- 
sical picture dSimon & Leighton| [l964) is a linear thermal con- 



vection scenario which associates supergranulation cells with 
the second ionisation of helium. Other scenarios put forward 
the idea that supergranulation is a surface phenomenon, as ac- 
tually indicated by th e recent results of loca l helioseismology 
dGizon & Birchl 120051) . IRieutord etail d2000h suggested that it 
could be the re sult of a large-scale i nstabi lity of surface granu- 
lar convection. iRincon & Rieutordl (l2003h proposed that some 
fixed flux boundary condition is imposed by the granulation 
to the layers just below: with these boundary conditions, the 
buoyancy destabilises the largest available scales. The associ- 
ated convective mo tions have a fairly low intensity contrast 
(Ivan der Borght 1974), very much like s uperg ranulation cells 
(Meunier et alJ2007l) . lRincon & Rieutordl (120031) further showed 
that an effective finite but large horizontal convection scale can 
be obtained in this framework by considering the dynamical ef- 
fects of a mean magnetic field on the flow. Collective plume in- 
teractio ns dRastl 120031: ICrouch et alj|2007l) and travelling- wave 
sheared dGreen & K osovichev 2006) convection or magnetocon- 
vection dGreen & Kosovichevl2007l) have also been suggested as 
a possible origin of supergranulation. 
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Overall the constraints on the large-scale dynamics of solar 
convection imposed by models and theories are still very loose. 
Guidance from observations is thus much desired. Particularly 
if a dynamical connection exists between granulation and su- 
pergranulation, some hints of such a scale interaction should 
be found in the dynamics of intermediate horizontal scales. 
Flows in the 3-10 Mm range are little known. This range of 
scales is traditional ly referred to as the m esogranulation range 
after the work of iNovember et al.l (119811) . but the very exis- 
tence of genuine enhan ced convectiye mo t ions at these scales 
has been much debated ([Straus et al.lll992t IStraus & Bonaccinil 
1997tlRieutord et al.l|2000HShine et al.ll2000l) . As pointed out in 
Rieutord et al.l (120001) . the main problem with mesogranulation 
was its way of detection, namely through the measurement of 
the horizontal divergence of the flows, which is highly sensitive 
to the way data are reduced. The very nature of the flows at these 
scales is nonetheless important to answer the question of the ori- 
gin of supergranulation. 

In order to make progress, it is clear that the large-scale side 
of the granulation spectral peak should be investigated in more 
detail. Sim ple models like the mi xing-length theory or plume 
dynamics dRieutord & Zahn|ll997l) do not predict any special 
spectral feature when the scale (horizontal or vertical) grows. 
Hydrodynamic nu merical simulations like the recent ones of 
IStein et al.l (E009) seem to go in the same direction and do 
not exhibi t any spectral feature reminiscent of supergranulation 
(see also Nord lund et al.l l2009). However, kinetic energy spec- 
tra derived from the radial velocity measured by SOHO/MDI 
(i.e. from dopplergrams) do show a rise of the spectral kinetic 
energy density at a s cale of 11 Mm and a peak at 36 Mm 
dHathawavet aLl lToOO). These results are clearly confirmed by 
the recent measurement of the horizontal surface flows at disc 
centre wit h the wide-fiel d high resolution camera CALAS at Pic- 
du-Midi dRieutord e t al. 2008). These observations are therefore 
clearly at odds with the most advanced numerical models of so- 
lar surface convection. 

In this paper we wish to provide accurate observational facts 
to guide modelling efforts and impose more constraints for the- 
ories. For this purpose we use new data sets collected by the 
Hinode mission and especially by the SOT instrument. Because 
of their high spatial resolution and their very low noise, these 
data are appropriate to either refine previous measurements or 
investigate new properties of the flows. In particular the subgran- 
ulation range can now be investigated as never before. We shall 
focus our attention on the power spectra of the flows because 
the information that can be obtained from such representations 
is usually the most robust aspect of these highly turbulent flows. 

The paper is organised as follows: after a brief presentation 
of the data and the techniques to process them (Sect. [2]), we detail 
the spectral properties of the measured velocity fields (Sect. [3} 
and the intensity field (Sect. [4]). A discussion and physical inter- 
pretation of the results follows (Sect. [5]). The main results and 
conclusions are summarised at the end of the paper (Sect. [6]). 

2. The data and the reduction procedure 

We used multi-wavelength data sets of the Solar Optical 
Telescope (SOT) on board the H INODE^ mission (e.g. 
Ilchimoto et al] 12004 ISuematsu et al] 120081) . The SOT has a 

1 The Hinode spacecraft launched in 2006 was designed and is now 
operated by JAXA (Japanese Space Authorities) in cooperation with 
NASA (National Aeronautics and Space Administration) and ESA 
(European Space Agency). 
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Fig. 1. Velocity field derived from the motion of granules dur- 
ing a 45 min window, interpolated by a Daubechies wavelet of 
4*953 km. Contours mark the horizontal divergence of the ve- 
locity. 



50 cm primary mirror with a spatial resolution about 0.2"at 
550 nm. 

2.1. Horizontal velocities 

A first set of images taken at disc centre was recorded con- 
tinuously from 29 August 10:17 UT until 31 August 10:19 
UT 2007, except for an interruption of seven minutes on 30 
August at 10:43 UT. We used the set taken with the broad- 
band filter imager (BFI) with a spectral width of 0.22 nm 
( Wedemev er-Bohm & Rouppevan derVoorll [2009) in the blue 
continuum at 450.45 nm. The mean time step between two 
successive frames is 50.1 sec. The field-of-view with BFI is 
1 1 1'.'6 x 1 1 1'.'6 with a pixel of 07 109 (1024 x 1024). After align- 
ment, the useful field-of-view was reduced to 104'.'6 x 104'.'6 or 
76.3 x 76.3 Mm 2 during 33h. The data were also k - u filtered 
with a phase velocity threshold of 6 km s _1 to keep only convec- 
tive motions. 

We determined the horizontal velocity field for various time 
averages from this set of data, us ing our algorithm of granule 
tracking CST dRieutord et al .1120071) . We recall that because of 
the nature of granules their motion is representative of the large- 
scale plasma flow only for length scales larger than 2.5 Mm 
and time scales longer than 30 min. These limits have bee n de- 
rived from numerical simulations by iRieutord et al .1 (1200 ll) and 
are confirmed by these data (see below). 

Other data sets were processed in the same way to compare 
several flow characteristics (see Table 1). 

2.2. Vertical velocities 

To determine the vertical component of the velocity we used a 
second data set from SOT/Hinode, also taken at disc centre and 
recorded on 4 September 2009 21:08 UT to 22:25 UT. These 
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Date 


Instrument 


FOV (Mm 2 ) 


Filter 


Time slot 


/i S G (Mm) 


£s G (km j /i 2 ) 






Wide-field data 








1996 May- Jul 


SOHO/MDI 


global 


Ni I, 676.8 nm 


62d 


36.4 


(200) 


2000 Apr 22 


TRACE 


225x250 


WL 


6x3 h 


35.6 


200 


2007 Apr 24 


TRACE 


238x275 


WL 


2x3.4 h 


29.8 


490 


2007 Mar 13 


CALAS 


290x216 


575 nm 


2x3 h 


36.4 


350 






Small-field data 








2007 Mar 8 


Hinode/SOT 


75.9x76.8 


G-band 


4h 


no peak 


220 


2007 Mar 15 


Hinode/SOT 


77.3x33.9 


G-band 


3x3h 


57/17-27 


240/350 


2007 Aug 29-30 


Hinode/SOT 


76.6x76.6 


blue 


llx3h 


28.7 


600 


2007 Oct 1 


Hinode/SOT 


2x(76x76) 


blue 


2x3. 8h 


27/31 


460/545 




Fig. 3. Left: kinetic energy spectra of the horizontal velocity field for various time windows. The vertical line marks the 2.3 Mm 
scale below which noise dominates. The black dots mark the scale at which the spectrum of a given time average disconnects from 
the 30-min one. Right: relationship between the length scale and the correlation time scale for the horizontal velocity field. Here 
N x = N y = 80 and N p = 256. 



data were taken with the narrowband filter imager (NFI) with a 
spectral width of 9 pm, near the Fel line at Aq =557.6 nm, sam- 
pling nine wavelengths in the line profile. This regular sampling 
allowed us to accurately measure the Doppler shift of the line 
every 28.5s. The field-of-view of this series is 82" x 82", which 
represents 60 x 60 Mm 2 . The pixel size is 0708. 

These data are completed by a third data set from 
SOT/Hinode, taken at disc centre and recorded on 4 September 
2009 7:35 UT to 10:17 UT with a regular cadence of one im- 
age every 57.5 s. These data were also obtained with the NFI 
but around the Fel line at 525.0208 nm. However, in that case 
the line profile was sampled only at two wavelengths, namely 
Aq - 1.8 pm and Ao + 10.8 pm, which did not allow us to de- 
termine the Doppler shift of the line in the same way as with 
the second data set. The radial velocity was instead estimated 
through the simple relationship 



v = k- + v , 

h + h 



where Vo and k are constants and and //, are the intensities 
in red and blue positions of the filter (+10.8 pm and -1.8 pm) 
respectively. The actual relationship between the radial velocity 
and the intensity difference for this particular line and positions 
of the filter is shown in Fig. [4] Note that an unavoidable pitfall 
with NFI is that the red and blue positions of the filter are in 
general asymmetric with respect to Aq. From the previous plot 
we see that this implies that large velocities cannot be measured 
correctly on one side of the line for these data (here on the blue 
side). The field-of-view of this series is 112" 6 x 112'.' 6 with a 
pixel of (X'16, corresponding to a physical field of 82 x 82 Mm 2 . 

The Doppler fields were then k - <±> filtered to remove fast 
horizontal motions of the patterns with a threshold of 6 km s . 

Unfortunately no flat-field image was available for the data 
sets taken at Aq =557.6 nm. We tried to circumvent this problem 
by creating a flat-field image for each image. For this purpose 
we averaged 50 images around the desired one and smoothed 
the result with a window of 16"xl6". The operation possibly 
altered the signal at scales larger than 12 Mm. Note finally that 
unlike the Aq =557.6 nm data sets, we were able to correct the 
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Fig. 2. Mean spectra and fluctuations: eleven independent hor- 
izontal kinetic energy spectra (dotted lines) computed with a 
three-hours time window, along with their average (thick solid 
line) are represented. Note that the fluctuations may reach an 
amplitude of 40% compared to the mean value. The straight line 
on the right indicates the power law E(k) ~ k characteristic of 
decorrelated noise. The vertical straight line marks the 2.3 Mm 
scale. Here N x = N y = 80 and N p = 1024. 




Fig. 4. Relationship between radial velocities and relative in- 
tensity fluctuations, using the two wavelengths at -1.8 pm and 
+ 10.8 pm aside the Fel line at 525 nm. Positive velocities are 
towards the observer. 



data set taken around Aq =525.0 nm with a genuine flat-field 

image. 

Finally, following iDanilovic et alJ d2008l) we computed the 
PSF and MTF of SOT to deconvolve all the images including 
those used to derive the Doppler shift from the instrument effects 
(central obscuration, spider and ccd). This correction turned out 
to be crucial for the investigation of subgranulation-scale tur- 
bulence, which is one of the novelties brought by the Hinode 
mission. 




l/\ (Mm" 



Fig. 5. Kinetic energy spectra of the horizontal velocity field dur- 
ing the emergence of two pores in the field-of-view. The solid 
line corresponds to the first series, at the end of which one pore 
has emerged. The dotted one corresponds to the second series, 
during which the second pore emerges. For each series, two 
time windows were used. The vertical dotted lines mark the 
12.5-75 Mm range where the supergranulation peak is found 
(cf. Fig.E). Here N x = U,N y = 82 and JV P = 1024. 



3. Velocity fields 

3. 1 . Horizontal velocity fields 

We give in Fig. \T\ a typical view of the velocity field using a 
time average of 45 min, on which the supergranulation pattern is 
clearly visible. As is customary with any turbulent flow, it is cru- 
cial to investigate the spectral signature of these flows. The high 
quality of the data enables us to use the highest spatial and tem- 
poral resolution s available to granule tracking. Using the CST 
algorithm, as in lRieutord et al. d2008l) . we computed the kinetic 
energy spectrum Ei t (k) for various time windows for the hori- 
zontal motions as shown in Fig. [2] We recall its definition 



Eh{k)dk , 



where Vf, denotes the horizontal velocity field averaged over the 
time window during which the granule tracking is performed; 
should be an ensemble average; here it is a spatial average, 
assuming that the turbulent flow is ergodic. 

For comparisons with other data sets, we plotted the dimen- 
sional values (in km 3 /s 2 ) of the kinetic energy spectral densities 
in all the figures. Since this value depends on the normalisation 
of the Fourier transform, we give the details of the expressions 
used in the paper in Appendix[A] In order to give a better spectral 
resolution, all spectra were computed using zero-padding. Thus, 
we give the size of the data (N x ,N y ) and the size of the padding 
area N p for each spectrum. 

These spectra are interesting for various reasons. First, we 
notice in Fig. [2] the prominent supergranulation peak centred 
at 30 Mm. This length scale is s lightly shorter than the one 
measured by Rieutord et al. (2008), who found the maximum at 
36.4 Mm. This change in the peak wavelength does not seem 
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to be an artefac t of the smalle r field. Indeed, reducing the field 
on the data of Rieutord et al. (2008) to a size comparable to 
that of the present data, one still obtains a value of 36.4 Mm. 
The difference seems to result from an intrinsic variability. To 
verify this point, we analysed other data from various sources 
(other Hinode fields, some data from the TRACE satellite). The 
results, together with t hose obtained from SOHO/MDI data by 
lHathawav et ail d2000Pl are summarised in Table Q] These data 
indeed show some variability both in the location of the peak and 
its amplitude. The fact that the small fields-of-view tend to give 
smaller values of the supergranulation scale may not be statisti- 
cally significant. 

Figure [2] also displays a plot of the individual spectra to- 
gether with their average. These spectra provide further evi- 
dence of the strong variability of the horizontal velocity field 
even at the supergranulation scale. They are constructed from 
three-hours-averaged velocity fields. We note that they do not 
show any energy excess within the mesogranulation range, i.e. 
between 3 and 10 Mm. The absence of any specific spectral fea- 
ture in this range is confirmed when other time windows are used 
(Fig.0. 

In Fig. [3] (left) we plot the mean spectra obtained with time 
windows ranging from 0.5 to 24 hours. Each spectrum is an av- 
erage of all the independent spectra that could be extracted from 
the whole series. For instance the 0.5h time window allowed 
us to extract 66 independent spectra from the 33 hrs of data. 
Note that all spectra show a distinct rise at scales below 2.3 Mm 
(l/A > 0.43 Mm- 1 ). This rise follows a linear law (E h (k) ~ k) 
as expected from a two-dimensional uncorrelated noise. Hence 
we clearly find the spatial resolution limit imposed by granule 
tracking. This confirms the analysi s performed on the numerical 
simulations by Rieu tord et alJ (1200 ll) . which suggested that gran- 
ule tracking could not be used at scales smaller than 2.5 Mm. 

The 30-minutes spectrum is the closest to an instantaneous 
one. This spectrum displays three features: (i) the supergran- 
ulation peak, (ii) a minimum in the kinetic energy density at 
10 Mm and (iii) a slow monotonic inc rease towards th e small 
scales. Feature (ii) was also found by lHathawav et alJ £2000) 
in the Doppler measurements from SOHO/MDI; their spherical 
harmonic spectrum shows a minimum near I = 400 which cor- 
responds to A = 11 Mm, quite close to our value. 

We note that the spectra obtained with time windows ranging 
from 30 min to 24 h almost superpose at large scales. However, 
shifting the attention to smaller scales, one finds the existence of 
a scale A, where the spectrum disconnects from the 30 min one 
(marked with black dots in Fig. [3}. The time average acts as a 
low-pass filter which removes (or simply weakens) the flows on 
shorter time scales. The size of the time window can hence be 
used as a proxy of the correlation time of the flow at the scale 
marked by the dot. Figure [3] (right) shows the scale A, as a func- 
tion of the time average associated with the time window. We 
find that the relationship between this length scale and the win- 
dow size is close to a power law A, ~ t 2 ^ . As the correlation 
time is ~ A/v^, we can relate v,\ and A and find ~ A" 3 ^ 2 . This 
translates into a power law spectrum E(k) ~ k 2 for the kinetic 
energy. This exponent is typical of the range 7-25 Mm. 

Finally, while processing other data sets to better evaluate 
the variability of the supergranulation peak, we came across the 
data recorded in the G-band by Hinode/SOT on 8 March 2007. 
Two series of images taken at [2:55, 5:45] UT and [6:00, 9:52] 



2 The amplitude of the peak for these data, given in between paren- 
thesis, is an estimate resulting from a crude conversion of Doppler sig- 
nal to horizontal velocities. 



UT drew our attention. As shown in Fig. [5] the supergranulation 
peak has disappeared if we consider the shortest time sampling. 
In this figure we distinguished the two series because one pore is 
emerging at the end of the first series (around 4:40 UT) and the 
second pore emerges in the second series (around 6:18 UT). We 
note that the kinetic energy density in the supergranulation range 
has weakened at scales shorter than 33 Mm in the second series 
when the pair of pores has reached its steady state. We also ob- 
serve that the TRACE data taken in April 2000 during the solar 
maximum also show a weaker supergranulation peak compared 
to the analogous data of 2007 taken at the solar minimum. Hence 
we note a trend that magnetic fields seem to markedly affect the 
amplitude of the supergranulation flows. This is likely connected 
with the anti-correlation between the size of the supergranules 
and the stren gth of the magnetic field observed by Me unier et al.l 
(120071120081) . 

3.2. Vertical velocity fields 

Let us now focus on the vertical velocity fields v z (x,y). As men- 
tioned above, these are obtained from the Doppler effect and thus 
each pixel of the image gives a measure of this component of the 
velocity field. We now compute the horizontal spectral density of 
the vertical velocity, namely 

i r°° 

-(v*) = E v (k)dk 

for various time averages. 

Figure|6]shows the spectra of the Doppler radial velocity ob- 
tained with the bisector method after reconstruction of the line 
profile at 557.6 nm. We extracted velocities at two different alti- 
tudes. A line chord of 16 pm was used to sample a layer a round 
an altitude of 80 km dAltrock et alJl975tlBerrilli e t al. 2002]) and 
a second chord of 8 pm, deeper in the line, was used to sample 
the atmosphere around 190 km. In order to detect potential spu- 
rious effects associated with the artificial flat-fielding operation, 
an additional set of spectra was computed from images not flat- 
fielded (Fig. 13). 

First these spectra seem to follow E v (k) ~ k 2 in the 2.5- 
7 Mm range. If the flat-fielding operation is correct, this power 
law may extend up to 30 Mm. Note that this slope is reminiscent 
of the power law found from the relationship between correla- 
tion time and the scale of horizontal velocities (see the previous 
section and Fig.[3]-right). In terms of velocity, the k 2 power law 
for the spectral density implies v. ~ y/kE v (k) ~ k 3 ^ 2 ~ A~ 3 ' 2 . 

As noted above, the quality of the data allowed us to inves- 
tigate the subgranulation dynamics accurately for the first time. 
The absence of Earth atmospheric turbulence and the deconvolu- 
tion of the images by the instrumental transfer function leads to 
high quality horizontal spectra of the vertical velocity. As shown 
by Fig. |6l we observe in this range of scales a decrease of the 
spectral density closely following a Ar 10 ^ 3 power law. When the 
signal is slightly averaged over a 4-min interval, the cut-off is 
steeper and the exponent close to -4. 

Concerning the granulation peak, the spectra show that the 
privileged scale is 1 .7 Mm. The spectral kinetic energy density 
is around 100 km 3 /s 2 at a mean altitude of ~80 km. At an al- 
titude around 190 km, the same spectra show an amplitude of 
50 km 3 /s 2 . This reduction of the spectral kinetic energy density 
is likely associated with the braking of upward motions which 
hit the stably stratified atmosphere. 

To confirm these results with independent data, we used the 
third set presented in Sect. 12.21 taken with NFI around the Fel 
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X (in Mm) : 30 15 10 7 5 



0.5 



100.00 



10.00 



1.00 



0.10 



0.01 






2 min 
4 min 
8 min 
14 min 
36 min 
72 min 



X = 557.6 nm 
Date = 4/9/2009 



alt. 80 km 

I 



0.01 



0.10 
1/X (Mm" 1 ) 



1.00 



X (in Mm) : 30 
100.00 



10.00 



15 10 7 5 



0.5 



1.00 



0.10 



0.01 
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8 min 
14 min 
36 min 
72 min 



X = 557.6 nm 
Date = 4/9/2009 



alt. 190 km 



0.01 



0.10 

1/X (Mm" 1 ) 



1.00 



Fig. 6. Kinetic energy spectra of the vertical velocity field for various time windows. The vertical dotted line marks the maximum 
of the spectrum at A — 1.7 Mm. Left: spectra computed with Doppler shifts taken with a line chord of 16 pm, corresponding to an 
mean altitude of 80 km. Right: same as left but for a chord of 8 pm, corresponding to an mean altitude of 190 km. The original 
images were flat-fielded according to the procedure described in Sect. 12.21 Here N x = N v = N p = 1024. 



line at 525.0 nm. Unlike the previous data at 557.6 nm we only 
had access to two wavelengths around the Fel line and had to 
derive the Doppler shift by the poor man's method, following 
Eq.Q] The results, shown in Fig. [8] are still interesting. With this 
method a k 2 power law is also obtained at large scales both on 
the 557.6 nm and 525.0 nm data. The latter data, which benefit 
from a good flat-field image, also indicate that the k 2 power law 
extends up to 30 Mm. 

However, we find that the cut-off is sharper on the subgran- 
ulation side, not far from a k^ 11 ^ power law. This is presumably 
a residual of the small-scale cut-off of the intensity power spec- 
trum, which follows the same power law (see below). This result 
suggests that the poor man's method is probably not appropriate 
for the computation of the small-scale radial velocity field. 

4. Intensity spectra 

With the data sets presented in Sect. 12.11 and 12.21 we finally in- 
vestigated the power spectrum of the horizontal intensity fluctu- 
ations defined according to 

1 / 2\ f 00 

-p| = J E,(k)dk. 

The spectra obtained from the blue continuum images at 
A = 450.4 nm are shown in Fig. [9] The variations in the 
mesoscale range are characterised by a spectrum close to a k 2 
power law, remarkably similar to the spectrum of vertical veloc- 
ities. Assuming that the surface radiates like a blackbody, these 
fluctuations reveal those of temperature. 

At scales smaller than a maximum located around the scale 
of 1.5 Mm, the power spectrum decreases steeply. In the sub- 
granulation range, the spectral density almost follows a power 
law in k~ 3 . We note that these results a g ree w ith those of 
IWedemever-Bohm & Rouppe van der Voortl (12009) as derived 
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Fig. 7. Same as in Fig. [6]- left, but with data not flat-fielded. 



from their Fig. 9. Our longer time series provides a better statis- 
tic, which probably improves the determination of the large- 
scale side of the spectra. 

We then used the images taken at 525.03 16 nm with the NFI, 
whose width is 9 pm. It turned out that this light come s from the 
contin uum just besides the Fel line at 525.0208 nm (Uitenbroek 
2009). As shown in Fig. [9] (right), the intensity spectra show a 
steeper cut-off in the subgranulation range than for the previous 
data, close to a k~ ll/3 law. We note that the granulation scale 
peak is at 2 Mm for this wavelength. 
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Fig. 8. Left: same as in Fig. [7] but with Doppler velocities computed with just two wavelengths in the line. Right: same as left but 
for A = 525.0 nm (images are corrected using a flat-field image provided by the Hinode mission, N x = N y = 704 and N p = 1024). 




Fig. 9. Spectral energy density of the intensity field. Left: data from the BFI at 450.4 nm (spectral width of 0.22 nm sampling 
the blue continuum) for various time windows (data set of 29-31/8/2007). The spectra are normalised so that the amplitude of the 
granulation peak is unity. The vertical dotted line marks the position of the granulation peak determined with instantaneous spectra. 
Here N x = 941, N y = 934 and N p = 1024. Right: data from the NFI at 525.0208+0.0108 nm (spectral width of 9 pm sampling the 
continuum on the red side of the line); here N x = N y = 704 and N p = 1024. 



Finally, the images taken at the blue side of the 525.0 nm Fel 
line are in fact very clos e to the line core, na mely at 525.019 nm. 
Atmospheric models of lUitenbroekl J2009) suggest that this in- 
tensity emerges at an altitude of 200 km above the continuum. 
The spectra shown in Fig. [TUl show a clear shift of the granula- 
tion peak towards larger scales, namely 2.5 Mm. No power law 



arises in the meso-scales, while the subgranulation range shows 
a roundish cut-off with no clear power law. 
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Fig. 10. Same as Fig. [9] -right, but with the NFI almost centred 
on the line core at A — 521.019 nm. 



5. Discussion and physical interpretation 

5.1. Large scales 

The foregoing results first confirm that the power spectra of hor- 
izontal flows at the Sun's surface are dominated by two scales 
corresponding to granulation at ~1.7 Mm and superg ranulation 
at ~30 Mm. This result confirms the one of lHathawav et alj 
(2000j), who used only one velocity component (the radial ve- 
locity) of the flow field. As far as the vertical velocity is con- 
cerned, the granulation scale clearly emerges at 1.7 Mm. On the 
large-scale side of granulation, we found that the vertical veloc- 
ity field decreases like A~ 3 ^ 2 with the horizontal length scale A, 
corresponding to a £ 2 -law for the associated kinetic energy spec- 
tral density. This law seems to be valid up to the supergranula- 
tion scale since its prediction of a supergranule vertical velocity 
corres ponds to the value of 30 m s _1 found by Hath away et alJ 
(2002) with MDI/SOHO data. Note that this conclusion holds 
only if the flat-field operations on the images did not damage the 
data. 

The £ 2 -scaling law of vertical kinetic energy can be related to 
the scaling law of the horizontal kinetic energy spectral density 
using physical arguments. Indeed, mass conservation (see below 
Eq. [2J implies that v z /H ~ v/JA, where H is the vertical scale 
height. If H is independent of A, and since v z ~ /T^ 2 , then Vh ~ 
/T 1 / 2 , meaning that E/,{k) ~ const. This law is approximately 
consistent with the horizontal flow data: as shown by the 30- 
minutes spectrum in Fig. [3] the spectral density Eh(k) does not 
vary very much in the 2.5-20 Mm range, where it stays between 
190and260km 3 /s 2 . 

To conclude on the k 2 power law, we note that it also appears 
in the intensity fluctuation spectrum obtained with the broadband 
filter in the blue continuum. This looks like a strong correlation 
between buoyancy fluctuations traced by intensity fluctuations 
and vertical velocities. However, we observe that this scaling is 
not as conspicuous at another wavelength with a narrower filter, 
which likely samples different layers (see discussion below). 

The derived horizontal flows may also be used to obtain 
some information on the scale height of the variations of the 
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Fig. 11. Schematic view of the results. We indicate the layers 
where power laws are perceptible in the subgranulation range. 



vertical velocity at the scale of supergranulation. Assuming that 
the anelastic approximation can be used at these scales (it means 
that acoustic waves are filtered out), we can write the mass con- 
servation as 



dz v z - ~ v zd z In p - V/, • v/, 



(2) 



where Vf, is the horizontal velocity fie ld and p is th e backg round 
density. Using the simple model of Simon et a 

D (fl99lf) for a 

horizontal divergent source, namely v r = vr/Rexp{-r 2 /R 2 ), 
we could estimate for a well-formed superg ranule that the di- 
vergen ce is approximately 10~ 4 s _1 . From Istein & Nordlundl 
(1 19981) . d z lnp ~ 2.4 MirT 1 near the surface. Then, using the 
estimate of 30 m s _1 for the vertical velocity at this scale, we 
find that the vertical scale height is ~ 1 Mm, which indicates 
that the supergranulation flow is likely shallow. This qualitative 
estimate agrees wit h recent local h elioseismic inferences based 
on Hinode data dSekii et al J 12007). which indicate that the im- 
print of supergranulation flows disappears around 3 Mm below 
the surface. Note that the foregoing exercise has been done us- 
ing typical values of all quantities and just gives a first impres- 
sion. To go further and get a more reliable view, this calculation 
should be applied to a large sample of supergranules. 

5.2. Sub-granulation scales: an injection-diffusive range ? 

On the other side of the granulation peak, in the subgranulation 
range, the power spectrum of vertical velocity closely follows a 
£-10/3 j aw ^ w hj cn steepens to k~ 4 with time averaging, while the 
intensity variance spectrum seems to follow either a k~ 3 law or 
a k~ 11 ^ 3 law depending on the data set used. These results are 
summarised in Fig. [TT] 

A natural question to ask is what these scaling laws tell us 
about the physics of the solar surface. It is indeed well estab- 
lished that a variety of characteristic power law regimes can be 
derived for all kinds of turbulent complex flows depending on 
the detaile d physical processes at work in the flow (e.g. shear 
turbulence (lLumlevlll967LlLo"hse & Miiller-Groe ling 1996 ), sta- 
bly an d unst ably st ratified turbulence (Obukhov 1959; Bolgiano 
1 1 9591 1 1 962t |L' vovll 1 99 lh . MHD turbulence (ISchekochihin et all 



2008, and references therein)) and on the relevant range of 
scales (e.g. injection, inertial or dissipation range) for each phys- 
ical field in the problem. In the paragraphs below, we therefore 
attempt to determine if the observed subgranulation scaling ex- 
ponents betray a particular flow regime. 

First of all, considering the many differences between homo- 
geneous, isotropic, incompressibl e turbulence and solar surface 
turbu l ence at observ able scales (Nordl und et al.llT997t iPetr ovav 
l200lt lRjnconll2007h . there is no real reason to expect that the 
Kolmogorov phenomenology of turbulence (which predicts a 
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AT 5 / 3 power law for the velocity field in the inertial range) ap- 
plies in the vicinity of the granulation scale. 

A possible interpretation for the joint observation of the 
-10/3 and -17/3 exponents is in terms of buoyant dynamics 
in the injection range of solar surface turbulence (i.e. not in the 
inertial range). Let us assume that (i) buoyancy fluctuations are 
directly proportional to temperature fluctuations (i.e. pressure 
fluctuations are not influential) and that (ii) the dynamic is dom- 
inated by a balance between inertial terms (or acceleration) and 
buoyancy forcing in the momentum equation. Further assuming 
that velocity and temperature perturbations have the following 
horizontal scale A dependence 

v,,(A) ~ A a , v,(A) ~ AP, 6T(A) ~ A y 

and that vertical variations are controlled by the typical pressure 
or density scale height H independent of A, mass conservation 
imposes that v z /H ~ Vf,/A, implying a = (5+1. Assuming finally 
the locality of non-linear interactions, (i) and (ii) translate into 
the dynamical balance 

v h v z /A + v 2 JH ~ ST , 

hence y = 2/3. Now, if the horizontal fluctuations of tempera- 
ture follow a Ar 17 ^ 3 spectrum, i.e. ST ~ V/t 17/,3_1 ~ A 1 ^, then 
y = 7/3, /3 = 7/6 and the horizontal spectrum of the vertical 
velocity is E v (k) ~ Ar 10 ^ 3 , just as observecfl Hence it is possible 
to relate the horizontal fluctuations of intensity at the continuum 
level to those of the vertical velocity in the upper layers (80 km 
or 190 km above) in the subgranulation range. 

Still, this is only one half of the story, since the argument 
does not predict the -17/3 slope for the intensity spectrum (it 
uses it as an input). This is because two dominant balances are 
required to determine scaling laws consistently in active scalar 
turbulence. In the context of thermal convection, the second bal- 
ance is obtained by isolating the dominant physical effects in 
the tempe rature equatio n or, to be more precise, in the Yaglom 
equation dRinco n 2006), which is the analog of the von Karman- 
Howarth equation for the statistics of scalars. In the Boussinesq 
approximation (which is not suitable for the problem of solar 
convection), a kr 11 ^ spectrum for temperature fluctuations can 
only be obtained if (i) the following low Peclet number dominant 
balanc^3 

v z d z {T) ~ kA6T (3) 

is assumed in the temperature equation (where k is the ther- 
mal diffusivity), i.e. temperature fluctuations are the diffusive 
response to velocity fluctuations advecting a mean, slowly vary- 
ing temperature profile (T) and (ii) the velocity fluctuations 
are not influenced by temperature fluctuations and have a k~ 5/3 
spectrum (Kolmogorov inertial range), i.e. v. ~ A 1 ^ . Then, 
ST ~ A 2 v z ~ AJ' 3 , or equivalently Ej(k) ~ k~ 11 ^ . This is the 
so-called inertial-diffusive regime (see e.g. lLesieurl [r990). Note 
that for any velocity spectrum slope steeper than -5/3 (say -10/3) 
Eq. (0 predicts a spectrum much steeper than kr 11 ^ for tem- 
perature fluctuations. Hence, it remains to be understood at this 

3 Assuming isotropy, the dynamical balance would instead read 
v\lA ~ ST, leading to E v (k) ~ Ar 13/3 for E,(k) ~ Ar 17/3 .This argu- 
ment can be formalised using a generalised exa ct von Karma n-Howarth- 
Kolmogorov equation for thermal convection (Rincon 2006). 

4 The Peclet number of a thermal eddy is the ratio between its thermal 
diffusion time and its typical turnover time. It is effectively of the order 
of one or smaller in the sub-granulation range. 



stage how the two slopes observed at the solar surface can be 
made fully consistent using dynamical arguments. 

Another interesting result regarding the sub-granulation 
range dynamics is that the intensity fluctuations taken with the 
BFI in the "blue continuum" at 450 nm have a much weaker 
spectral cut-off close to a AT 3 power law (AT 4 if slightly aver- 
aged) than those taken in the continuum near the Fel line at 
525 nm. This is quite surprising, because the fluid layer con- 
tributing most to the blue continuum should be quite similar to 
the one contributing to the 525 nm continuum from which the 
AT 17 / 3 spectrum was derived. Beside the simple wavelength dif- 
ference, which indicates a slighly lower contributing layer for the 
450 nm radiation, the most obvious difference between the two 
signals comes from the spectral width of the filters. The width 
of the NFI of 9 pm is such that the continuum near 525 nm is 
"pure", while the width of the BFI at 450.45 nm is 220 pm and 
includes a few absorption lines. For this reason, the layers con- 
tributing to the light collected in the two wavebands could be 
somewhat different. At the moment we have no obvious expla- 
nation of these two different cut-offs for the intensity spectrum in 
the continuum. Of course we cannot totally exclude some unex- 
pected bias in the data processing, but if the result is confirmed 
by independent approaches, detailed numerical simulations of 
the dynamics are likely the only way to elucidate this puzzle. 

For the sake of completeness, we mention that a AT 3 spec- 
trum could be a simple signature of a physically smooth temper- 
ature field, which can be expanded in Taylor series. Indeed, if 
T(A) = T(Q) + AV h T u=0 then ST = T(A) - T{0) ~ A, or equiva- 
lently Ej(k) ~ AT 3 . But, similarly to the argument presented pre- 
viously, we underline that this scaling is not fully consistent with 
a Ar 10 ^ 3 scaling for the velocity field if a dominant balance be- 
tween inertia and buoyancy is assumed in the momentum equa- 
tion. 

The previous theoretical arguments should of course not be 
taken at face value. They simply suggest that it might be possible 
to construct some sort of "injection-diffusive" turbulent regime 
characteristic of the sub-granulation range, and serve to illus- 
trate that the standard framework of incompressible homoge- 
neous isotropic turbulence is probably too restrictive to theoret- 
ically explain the observations of the thermal boundary layer at 
the solar surface. 

These new observational measurements for the subgranula- 
tion range definitely call for a deeper investigation. A detailed 
theoretical understanding of the solar power spectrum at the be- 
ginning of the subgranulation range 0.1-1 Mm might be ob- 
tained by re -analysing existing num erical simulations of solar 
granulation (IStein & Nordlundll 19981) . lookingfor the dominant 
terms in the temperature (or energy) equatioifl Another possi- 
bility would be to perform idealised high-resolution numerical 
simulatio ns of thermal turbulence done in the low-Peclet approx- 
imation (Thua3[1992t lLignieresl fi999) to better understand this 
complex and fairly unusual turbulence regime from the theoreti- 
cal point of view. 

6. Conclusions 

The power spectrum of surface flows at the Sun's surface 
has been investigated using mainly image series from the 
Hinode/SOT instrument. The results may be summarised as fol- 
lows. 

5 It is remarkable that some of the early numerical simulations by 
Nordl und et alj fl997) seemed to produce a A -17 ' 3 scaling for intensity 
fluctuations. 
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1. The smallest scales, i.e. the subgranulation range, can be in- 
vestigated for the first time without a seeing effect and using 
corrections of the MTF of the instrument. On the one hand, 
the spectral cut-off of intensity fluctuations is similar to a 
kr 11 ^ power law for images taken with a narrow filter (width 
~ 9 pm) in the continuum on the red side of the 525 nm Fel 
line and to k~ 3 for images taken with the broadband filter 
(width ~ 220 pm) around 450.45 nm. On the other hand, 
the vertical velocity horizontal spectrum measured 80 km or 
190 km above the continuum at 557.6 nm shows a steep cut- 
off in AT 10 / 3 , which can be related to the AT 17 / 3 power law by 
invoking a balance between buoyancy and inertia. The val- 
ues of the two scaling exponents suggest that the temperature 
fluctuations impact the dynamics at a scale well below the 
one determined by the unit Peclet number. Finally, we also 
note the decrease of the granulation peak amplitude with al- 
titude, possibly betraying the deceleration of the upflows. 

2. The granulation spectral peak at 1.7 Mm is conspicuous in 
the vertical velocity spectrum. The location of the inten- 
sity spectrum peak depends on the wavelength: it is close 
to 1.5 Mm in the blue continuum at 450 nm, 2 Mm in the 
continuum near 525 nm, and to 2.5 Mm in the core of the 
Fel line at 525.0 nm, which probes the dynamics in a layer 
200 km above the continuum. 

3. Moving to the large scales, we note that the spectral density 
of vertical kinetic energy decreases as k 2 with decreasing k 
(equivalently, the amplitude of v, decreases as /t~ 3 ^ 2 with the 
horizontal scale A). This k 2 power law also emerges from the 
relationship between the scale and the correlation time of 
the horizontal flows. Furthermore we note that the spectral 
density of intensity fluctuations in the blue continuum also 
scales like k 2 in the 2.5-10 Mm range. 

4. The power spectrum of horizontal flows has been explored 
in the 2.5-40 Mm range. It clearly reveals the supergranu- 
lation peak as well as the large-scale side of the granula- 
tion peak. These observations c onfirm in passing the result 
of the numerical simulations of iRieutord et al.l (1200 ll) that 
granule tracking cannot determine plasma flows at scales be- 
low ~2.5 Mm. The energy density minimum is found at a 
scale of 12 Mm. These results confirm with a completely 
di fferent technique t he sph erical harmonics spectra obtained 
by Hatha way et alJ (l2000l) from dopplergrams. Focusing on 
the supergranulation peak, we could estimate its amplitude 
to be in the range of 200-700 km 3 /s 2 and show its strong 
variability. From Figs. [3] and [5] we may add that the visible 
part of the granulation peak points to an amplitude in the 
range of 200^1-00 km 3 /s 2 . Interestingly enough, we discov- 
ered that the emergence of a pair of pores is able to erase the 
supergranulation peak completely, which points out the sen- 
sitivity of the supergranulation amplitude flow to the ambient 
magnetic field. Finally, combining all the data, we tentatively 
estimated the scale-height of the vertical velocity field at the 
supergranulation horizontal scale. We found a scale height 
of approximately 1 Mm, indicating that supergranules are 
likely surface flows. However, this estimate should be ver- 
ified for a number of supergranules to obtain a statistically 
robust value. 

These observational results provide a set of landmarks that 
should be recovered by future numerical simulations. Future ob- 
servational work will focus on the large-scale side of the super- 
granulation peak to determine what scaling law, if any, controls 
the very large-scale dynamics of the Sun between supergranula- 
tion and differential rotation. Another point deserving a detailed 



investigation is the dynamical influence of the magnetic fields 
on the surface flows and their power spectrum. 
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Appendix A: Absolute spectra 

In this paper, we give the kinetic energy density spectra with 
their absolute values in km 3 /s 2 . To facilitate future comparisons, 
we document the way these spectral densities are computed in 
detail. We recall that the kinetic energy spectrum E(k) is defined 
as 



E(k)dk. 



(A.l) 



It is related to the Fourier transform of the velocity components 
by 



(v*(k)v j (k')) = 4> ij (k)6(k-k'), 



(A.2) 



where cf>ij(k) is the Fourier transform of the two-point correlation 
function of the velocities, namely, Qij(r) = (vi(x)vj(x + r)^. 6 is 
the Dirac distribution. In two dimensions (our case here), E(k) 
and (pij are related by 



E(k) 



i, r 2n 



(k)de k , 



(A3) 



where 0* is the polar angle of k. The numerical 2D-Fourier trans- 
form (like in the IDL software) is defined as 



-2in{k x x+kyy) 



N x N y 



(A.4) 



which means that the Dirac distribution in a discrete way reads 



= 1 V 



-2i7i(k x x+k y y) 



(A.5) 



Dirac distribution. Using the preceeding normalisation for both 
v and 6{k), we find that 



lv| 2 |v| 2 (AyV v )V(27r) 2 , , , 

L J = \v\ 2 N x N y p 2 l{2nf . (A.9) 



6{k) (N x N y )pH2n) 
If we restrict our case to a square data set (N x = N y = N), then 

where we noticed that the dimensional wavenumber k and the 
non-dimensional one are related by a factor 2n/Np. Finally, 



E ^ = \tT\ X Mk')\ 2 + \v Y (k')\ 2 , 



(A. 10) 



We[k,k+dk]/27r 



where v x and v y are the discrete Fourier transform of v x and v y . 
Note that the summation zZvk'e[k,k+dky2n implies that the k' fac- 
tor of the surface element k'dOkdk' is aut omatic ally taken into 
account. Sometimes the quantity P(k) = -^kE(k) is used; here 



'ir,|2 



k'\v\ 



(A. 11) 



Vk'e[k,k+dk]/2n 



Effect of zero padding: the Fourier transforms are done on a 
square of size N x N, while the data usually occupy a rectangle 
of the size of N x xN y (N x < N,N y < N); the energy of the Fourier 
modes is thus a fraction N x N y /N 2 of the total. Hence, to obtain 
spectra which are independent of the zero padding, we need to 
multiply the resulting spectral densities by N 2 /(N x N y ) so that 



E(k) 



N 3 p 



4nN x N y w ^ kk+dk]l2n 



Z 



(A. 12) 



Effect of the normalisation of the Fourier Transform: if the 

Fourier Transform is multiplied by a factor A, then E(k) is mul- 
tiplied by A and P(k) is multiplied by a factor a/A. The P(k) 
spectrum has no absolute significance. 



Now, spectral quantities in turbulence theory usually use 



/(*) 



f(r)e-' kr dxdy 



(A.6) 



as the definition of the 2D-Fourier transform. Setting k = 2nk' 
then 



/(*) = (2tt)- 2 j f(r)e- 2M ' r dxdy, 



(A.7) 



The correspondence is accordingly 

f(k) = f(k')N x N y p 2 l(2n) 2 and 8(k) = 6N x N y p 2 l(2n) 2 , (A.8) 

where N x and N y are the size of the data, p is the linear size of the 
pixel. We wish to derive the correct normalisation of E(k) which 
is that of (pu(k). We note that <pu(k) is the ratio of the squared 
norm of the Fourier components of the velocity divided by the 



